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o : 

O ' Abstract 

<N ■ 

^ . We compute the spectrum of the trigonometric Sutherland spin model of BCn type 

| in the presence of a constant magnetic field. Using Polychronakos's freezing trick, we 

Q . derive an exact formula for the partition function of its associated Haldane-Shastry 

spin chain. 

Q\ ; 

1 Introduction 

oo ; 

In 1971, F. Calogero [0] introduced a solvable quantum model describing a system of 
N particles with two-body interactions depending on the inverse square of the particles' 
distance. In the same year, B. Sutherland |22j proposed a similar model with an interaction 
potential of trigonometric type. The importance of these models, which was already 
apparent from the very beginning, is now widely acknowledged by the theoretical and 
mathematical physics community. From a mathematical point of view, these models are 
integrable both in their classical and quantum versions, in the sense that they admit a 
complete set of integrals of motion. Moreover, they are also exactly solvable, in the sense 
that their spectrum and eigenfunctions can be expressed in closed form. From a more 
physical point of view, Calogero-Sutherland (CS) models play an important role in many 
different fields, like for instance Yang-Mills theories |17) . the quantum Hall effect pQ, 
5-j \ random matrices |2fi| I28j and fractional statistics |18| 121]. 

Calogero-Sutherland models were cast into a very elegant mathematical framework by 
Olshanetski and Perelomov in [2U]. They showed that these models are limiting cases of a 
more general one with a two-body interaction potential of elliptic type, and uncovered their 
relation to the An root system. In fact, these authors also constructed generalizations of 
the previous models associated with all the classical (extended) root systems, like BCn. 

The extension of CS models to particles with internal degrees (typically interpreted as 
spin) of freedom was first proposed in the classical case in Ref. [Ej- in the quantum case, 
spin CS models were actively studied in the last decade, both in the An and BCn cases. 
It turns out that these models inherit the basic properties of their scalar counterparts, 
namely their integrability and exact-solvability. There are essentially two approaches in 
the study of spin CS models, namely the supersymmetric formalism jSJ |H1 EH an d the 
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Dunkl or exchange operator method El EH HH H2 1221 EH] • Usually, but not always 
the interaction of the spins with a constant external magnetic field is not included in the 
spin CS Hamiltonian. 

In 1988, Haldane and Shastry [23] independently introduced a new type of solvable 
spin chain with long-range position-dependent interactions. The sites of this chain are 
equidistant points in a circle, and the interaction between the corresponding spins is 
proportional to the inverse square of their chord distance. Shortly afterwards, Fowler and 
Minahan ^3] showed that this chain is completely integrable by means of Polychronakos's 
exchange operator formalism [22]. The connection with the Sutherland model, although 
already noted by Shastry in his paper, was made precise by Polychronakos in [23], using 
what he called the "freezing trick". The main idea behind this method, which can actually 
be applied to any spin CS model, is to take the strong coupling constant limit in the 
Hamiltonian, so that the particles become "frozen" at the equilibrium positions of the 
scalar part of the potential. In this way one can obtain new spin chains of Haldane- 
Shastry (HS) type, in which the sites are not necessarily equally spaced. Most of the 
literature on HS spin chains is devoted to those based on dynamical spin models of An 
type, while their BCn counterparts have received comparatively less attention. Yamamoto 
and Tsuchiya proved the integrability of the rational HS chain of BCn type (HOI) although 
they did not compute its spectrum. The trigonometric BCn spin chain was discussed by 
Bernard, Pasquier and Serban [3] in the spin 1/2 ferromagnetic case, but only for equally 
spaced sites. The integrability of the trigonometric/hyperbolic version of this chain was 
established in [ZlElj) although again its spectrum was not computed. In a recent work 
jl()j . we have extensively studied the trigonometric BCn spin chains (both ferromagnetic 
and antiferromagnetic) for arbitrary spin and without assuming that the sites are equally 
spaced. The main result in the latter paper is the derivation of a closed-form expression for 
the partition function of the model, following a method based on Polychronakos's freezing 
trick. 

In this paper we study the effect of the presence of a constant external magnetic field 
in the trigonometric BCn Sutherland model for spin 1/2 and in its associated spin chain. 
This modifies the Hamiltonian of the spin chain by the addition of a term proportional 
to the projection of the total spin operator along the direction of the magnetic field. It 
turns out that both the dynamical spin model (if the magnetic field is suitably oriented) 
and its associated chain remain solvable, although the method used in ^01 to compute 
the spectrum of the dynamical model, based on expressing the Hamiltonian in terms of 
a commuting family of Dunkl operators, cannot be applied in this case. We shall see, 
however, that the basis of the Hilbert space constructed in the latter reference to solve the 
model in the absence of a magnetic field can be slightly modified so that the Hamiltonian 
of the dynamical model is still triangular. From the spectrum of the dynamical model we 
shall then evaluate in closed form the partition function of its associated spin chain using 
the freezing trick. 

The paper is organized as follows. In Section 2 we introduce the dynamical spin model 
and its corresponding Hilbert space. The spectrum of this model is computed in Sec- 
tion 3, by showing that the term due to the magnetic field is triangular in an appropriate 
modification of the basis constructed in ^01- In Section 4 we define the HS spin chain 
and compute its partition function using the freezing trick and the results of the previous 
section. 
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2 Preliminary definitions 

In this section we set up the notation used throughout the paper and define the Hamil- 
tonian of the Sutherland spin model of BCjsr type in the presence of a constant external 
magnetic field. We shall denote by S the Hilbert space corresponding to the internal 
degrees of freedom of N identical spin \ particles. Let 

B C s = {\s) = \ Sl ,...,s N ) | Sl = U}, (2.1) 

be a basis of S whose elements are simultaneous eigenstates of the ( component of the 
one-particle spin operators, OQ being an arbitrary direction. The corresponding spin 
permutation and reversal operators Sij and Si are defined by 

Sij \s\, . . . , Si, . . . , Sj, . . . , S]\[) \S\ , . . . , Sj , . . . , Sj, . . . , S]\j) , 
Si\si, . . . , Si, . . . , Sn) = \si, . . . , —Si, . . . , sn) ■ 

We will also use the customary notation Sij = SiSjSij. The operators Si and Sij generate 
a multiplicative group isomorphic to the Weyl group of Bn type. Similarly, the coordinate 
permutation and sign-reversal operators Kij and Ki are defined by 

(Kij f)(xi, . . . , Xi, . . . , Xj, . . . , Xjv) = f{x\, . . . ,Xj, . . . ,Xi, . . . , Xjv) , 

(Kif)(xi, ...,Xi,..., x N ) = f(x u -Xi, x N ) , 

and Kij = KiKjKij. The total permutation and sign reversal operators will be denoted 
by Hi = KiSi and 11^ = Kij Sij. The multiplicative group generated by Ki and K^ 
(respectively, IT and 11^) is also isomorphic to the ^Ar-type Weyl group. 

Let us define the antisymmetrizer with respect to the symmetric group generated by 
IT j as 

~ 1 m 

■ i=i 

Pi being an element of this group and sgn(Pj) its signature. Likewise, we will denote by 
A e = 2- N H(l + eU t ) 

i 

the symmetrization (e = 1) or antisymmetrization (e = —1) with respect to sign reversals. 
Here and in what follows all sums and products run from 1 to N, unless otherwise stated. 
We shall make use of the projection operator 

A t = AA t 

on states antisymmetric under particle permutations and with parity e under sign-reversals. 
The Hamiltonian of the spin Sutherland model of BCm type in a constant external mag- 



(2.2) 
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netic field 25 is given by 

H e = — ^2 + a ^2 [ sni_2 ( x * ~~ x j) ( a + Sij) + sin _2 (xj + Xj) (a + Si 
+ b sin _2 Xi (b — eSi) + b' cos~ 2 Xj (b' — eSi) 



(2.3) 



e 2 



where £ = (£* E^, E c ) is the total spin operator and (0£, Or], OQ is an arbitrary system 
of orthogonal axes. Here we have assumed that the real constants a, b, b' are greater than 
5, and we have denoted by g and e the particles' gyromagnetic ratio and electric charge, 
respectively. The last (diamagnetic) term in Eq. ()2.3|) has to be dropped to preserve the 
solvability of H e . In fact, this term vanishes identically if the magnetic field is parallel 
to the x axis. Since our main interest is to study the spin chain associated with the 
Hamiltonian H e , we have preferred in what follows to drop the diamagnetic term in order 
to keep the direction of the magnetic field in the spin chain arbitrary. 

The above Hamiltonian possesses inverse-square type singularities for Xi ± Xj = kn and 
Xi = kn/2, with k £ Z. Since the nature of these singularities prevents the particles from 
overtaking each other and from crossing the singular hyperplanes X{ = kir/2, the particles 
may be regarded as distinguishable, with configuration space 

C= |x= (x!,...,x N ) eR N < xx < ■■■ < x N < |} . (2.4) 
The Hilbert space of the system can thus be taken as TL = Lq(C) <8> S, with 

L 2 (C) = {/ G L\C) 3 lim" |x l ±x i -fcvrr a |/(x)|, 3 lim" N~ 6 |/(x)| , 

l Xi±Xj— >K7T Xi— *U 

3 lim" |^-7r/2|- 6 '|/(x)|; A: = 0,1, 1 < i± j < n\ . 

X{— *7r/2 > 

The conditions imposed in the definition of Lq(C) guarantee the finiteness of (ip,H e ip) 
for all tjj E TL- It can be shown that the operator H t : 7i — > TL is actually equivalent to 
any of its extensions to spaces of symmetric or antisymmetric functions (with respect to 
both permutations and sign reversals) in Lq(C) <8> S, where C is the A^-cube (— f , and 
Lq(C) is defined similarly to Lg(C). For technical reasons, it is convenient to consider 
that H t acts in the Hilbert space 

H 6 = A t (L 2 (C) ® S) (2.5) 

of states antisymmetric under permutations and with parity e under sign reversals. 



3 Spectrum of the dynamical model 



The Hamiltonian (|2.3|) with 33 = was thoroughly studied in ^U], and its spectrum 
was exactly computed. The calculation was based on the fact that when 33 = the 
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Hamiltonian (j2.3fl is the image of the operator 

H' Q = — ^2 + a ^2 \ s ^ XL ~ 2 (. x i ~ x j) ( a ~ K-ij) + sin _2 (x, + Xj) + (a — K., 

i ijtj 

+ b^2 skT 2 Xj (b- Ki) + b' ^2 cos~ 2 Xi (b' - Ki) 



(3.1) 



under the mapping (cf. 

Kij i ^ i-» eSi . (3.2) 

The operator H' Q was expressed as a sum of squares of an appropriate set of commut- 
ing Dunkl operators Jj, i = 1,...,N, which preserve a flag IZq C 1Z\ C ••• of finite- 
dimensional spaces of smooth functions, where 

T^k = (<K X ) ex P (2i rijXj) | rij = —k, —k + 1, . . . , k , j = 1, . . . , N) , (3-3) 

j 

and 

(/>(x) = | sin(xj — xj) sin(xj + Xj)\ a | sinxj| 6 | cosxj| b . (3-4) 
For each value of k, we constructed a basis Bk of 7£fc such that the matrix representing 

oo 

Hq in this basis is triangular. Since the closure of U IZk is Lq(C), this observation 

fc=o 

immediately yields the spectrum of H' . The spectrum of H e in the absence of magnetic 

oo 

field was then computed by suitably extending the basis |J Bk to a basis B^ of the Hilbert 

fc=o 

space TL e ; see |TU] for the details. 

The operators Sij and Si can be expressed in terms of the usual one-particle spin 
operators Xli = (Ef,E^,E?) as follows 

Sij = 2Sj • Sj + — , Si = 2S| . 

Hence, if the magnetic field in the Hamiltonian (|2.M|) is directed along the £ axis, the term 
33 • 5] can be expressed in terms of the spin reversal operators Si as ^^X^Si, where 
05 = 1 35 1 . In this case, the Hamiltonian H t may be written as 



H e = — ^2 9^ + a ^2 [ s * n 2 ( Xi ~ x ^ ( a s * n 2 ( x * x ^ ( a 

+ 6 sin" 2 ^ (b - eSi) + 6' ^ cos" 2 ^ (6' - efli) - y 05 . (3.5) 

i i i 

Although this Hamiltonian is the image of the operator 

H' e =— ^2 ®xi + a ^2 [ srn_2 ( x * ~~ x j) ( a ~ Kij) + sin~ 2 (xj + Xj) (a — Ki 

i i^j 

+ b ^2 sin -2 ^ (b-Ki) + b'J2 cos~ 2 Xl (b' - K { ) - ^ 05 J2 K i > 
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under the mapping (|3.2j) . the Dunkl operator techniques of Ref. cannot be directly 
applied to prove its integrability, since it is not clear how to express the operator H[ in 
terms of a commuting family of differential-difference operators. Moreover, even if the 
new term proportional to the magnetic field leaves invariant the spaces IZ^, the matrix 

oo 

of H' e in the basis (J B^ of Lq(C) introduced in Ref. ^U] need not be triangular. It 

k=0 

is more convenient, therefore, to work directly with the spin Hamiltonian (|3,5|) and its 
representation in a slight modification of the spin basis B^ used in the latter reference. 

Let us recall, to begin with, the construction of the basis Be- To this end, we need 
to introduce the following notation. Let Nq denote the set of nonincreasing multiindices 
n = (n±, . . . with m = 0, 1, . . . and m > • • • > n^. If n,n' G Ng^, we shall say that 

n -< n' if ni — n'i = • • • = m—\ — n'i_i = and rn < n-. Let f n denote the function 

/ n (x) = 0(x)e 2i ^^, n£<, 

where ^(x) is given by (|3.4j) . The basis Be can be taken as any linearly independent subset 
of the set 

{A e (/ n |a» |nGN^, \s)eB C s } (3.6) 

ordered so that A e (/ n |s}) precedes A e (/ n /|s')) if n -< n'. It should be noted that the 
basis Bg in Eq. (|3.6|) could actually be replaced by any basis of the spin space S without 
changing the triangularity of the Hamiltonian H e \^_ . As shown in Ref. jlUj . a linearly 
independent subset of the set (|3.6|) is obtained by imposing the following conditions on 
the multiindex n and the spin basis element \s): 

1. = cardji | ni = m} < 2 — 5 m o for all m = 0, 1, . . . ; 

2. If m = TH+l, then Si = | and s i+ \ = |; 

3. If njy = 0, then sjy = f. 

It is more convenient for our purposes to work a slight modification B^ of the basis Be, 
obtained by replacing Bg in the previous construction by the basis 

B^ s = {\a) = \a 1 ,...,a N ) | ^ = ±1/2 } (3.7) 

of simultaneous eigenstates of the one-particle spin operators E|. Conditions 2 and 3 
above should accordingly be replaced by 

2'. If rii = rtj+i, then cij = +1/2 and 0*4.1 = —1/2; 

3'. If njy = 0, then = e/2. 

The last condition is due to the fact that the one-particle spin states 

|±l/2) = -L(|T)±|!» 
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have parity e under flipping of the ( component of the spin. Since the operator Yl% Si 
commutes with the projector A e , it is diagonal in the basis B>\. Indeed, 

J2s^A e (f n \a)) =A e J2fn(Si\a}) = (^2 f r i )A e (/ n |a)) = A(n, a)A e (/») . 

i i i 

Taking into account the conditions 1, 2', and 3' for the modified basis Be, it immediately 
follows that the eigenvalue A(n, a) can be expressed as 

A(n, a) = d + (n, a) — d_(rt, a) , 

where 

d±(n,a) = cardjz | #(rij) = 1 and Oi = ±1/2} . 

By the previous remark, the Hamiltonian H e \ <s _ is still triangular in the modified basis 
Bi, with diagonal elements 

El = J2 i 2n i + b + b' + 2a{N -i)f , (3.8) 

i 

cf. Ref [TJ3- It follows that the complete Hamiltonian (|3.5|) is triangular in the basis B\, 
with eigenvalues E na given by 

E na = E n -^j-<B\(n,a). (3.9) 

This formula for the spectrum of the Hamiltonian (J3.5|) will be used in next section to 
compute the partition functions of the corresponding spin chains. 

4 An HS spin chain of BCn type in a magnetic field 



Using Polychronakos's freezing trick |23| I24j . one may obtain a Haldane-Shastry spin 
chain of BCn type associated with the Hamiltonian ()3.5|) . This technique, thoroughly 
discussed in JOJ) consists in taking the large coupling constant limit a — > +oo, while 
maintaining constant the ratios (5 = b/a, f3' = b'/a and B = —eg^B/^Qa). In this limit, 
the eigenfunctions of H t become sharply peaked around a minimum of the potential 

U(x) = ^2 ( sin~ 2 (xi -Xj) + sin~ 2 (xj + xj)) + ^(/3 2 sin" 2 ar» + /3' 2 cos" 2 x») , (4.1) 

i^j i 

and thus the spin and the dynamical degrees of freedom decouple. It is important to note 
that there is unique minimum x° = {x\, . . . , x N ) of the potential U in the Weyl chamber 
C, as proved in Ref. [TO] . 
Let 



H s = - <9 2 . + a(a — 1) y] (sin 2 (xi - Xj) + sin 2 (xi + Xj)) 
i i^j 

+ 6(6 - 1) sin" 2 Xi + b'(b' - 1) ^ cos 



(4.2) 

2 
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denote the Hamiltonian of the scalar BCn Sutherland model. The Hamiltonian h e of the 
HS spin chain associated with H e is defined by 



h e = - (H e 
a 



namely 



sin- 2 (x° - a$ (1 + Sij) + sin^(^ + xfj (1 + Sy 



(4.3) 



2^0 



+ (P sin ~ 2 x i + 0' cos ~ 2 x i) (! - e <%) + 8B^2Si. (4.4) 



This Hamiltonian differs from the one in by the last term, which represents the 
interaction of the spins with a magnetic field along the £ axis of constant magnitude 
— 16B/(gq). As shown in Ref. Eq. (|4..'Sj) and the previous considerations lead to the 
relation 



Z e (T) = lim 



Z e (aT) 



a^oo Z s (aT) 



(4.5) 



between the partition functions Z e , Z t , and Z s of the respective Hamiltonians h e , H t , and 
H s . 

We shall now compute the partition function of the spin chain (|4.4|) in closed form by 
evaluating the RHS of Eq. (jO|) . From Eqs. (|3~H)) and (|33|) it follows that 



E ~ 



a 2 £ + 8a 



[/2 n i(P + N ~ i ) +BX{n,a) 



(4.6) 



where (3 = ^(/3 + /?'), and we have dropped the term independent of a which becomes 
negligible in the limit a — » oo. The constant Eq = 4 ^(/3 + A r — z) 2 , which is independent 
of n and a, can also be dropped from both partition functions Z e and Z s without modifying 
the value of Z e . With this convention, it was proved in ^01 that when a — > oo the partition 
function of the scalar Sutherland model ()4.2j) can be written as 



Z s (aT) ~ [] 



i(/3+JV-|(i+l)) 



-1 



(4.7) 



where we have set 

. = „-8/fc B T 



The calculation of Z e (aT) is more involved. To perform it, it is convenient to represent 



the multiindex n £ Nq appearing in (J4.6|) as 



n 



k± k 2 

mi, . . . ,TOi,m 2 , • • • ,m 2 , 



(4.8) 



where mi > • • • > m r > and /cj = #(?n.j) satisfies X/I=i ^» = N ■ The characterization 
of in the previous section implies that fcj G {1,2} and k r = 1 if m r = 0. We shall 
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denote by Vn the set of partitions k = {k\, . . . , k r ) of the integer N such that k% S {1, 2}. 
If k = (ki, ■ ■ ■ , k r ) S Vn, we define 

d(k) = card {i \k% = l} . 
Since d(k) = d+(n,a) + d_(n,er), it follows that 

X(n,a) = 2d + (n,a) - d(k) , (4.9) 
and therefore 

E na ~8a\j2MP + N -i)+ B [ 2 d+(n,<r)-d(k))] , (4.10) 

i 

where we have used Eq. Q4.6|) without the inessential ground state energy Eq. After 
expressing the first sum in terms of m and k as in Ref. ^U] we obtain the expression 



r 

E ~ 8a ' ' 



mi vi (k) + B(2d+(n,a) - d{k)) 



(4.11) 



i=l 

where 

The partition function Z e (aT) is therefore given by 

r 

Z e (aT) E q £™ iUiik) J2<l B[2d+{n ' (T) ~ m] > ( 4 - 12 ) 

kev N m ^>-> m r>o W) 

where the sum over the spins is restricted to those values of \a) such that A(f n \a}) 6 £>f. 
The latter sum depends essentially on whether m r > or m r = 0. Indeed: 

Case 1: m r > 0. 

In this case, for a given partition k € Vn, d+(n,a) can take any value in the range 
0, . . . , d(k). The number of states of the basis £>f for which d+(n, a) = 5 is given by the 
combinatorial number { d ^)- Hence 

M /^x . d{k) 



q B[2d+(n,a)-d{k)} = g ~Bd(k) ^ f d ^\ q ^BS = g -Bd(k) A + g 2B^ 

|<j> 5=0 ^ ^ ' 

\ <f(fc) 

g B + g" B , (m r >0). (4.13) 



Case 2: m r = 0. 

Note, first of all, that Condition 1 on the basis £>f implies that in this case k r = 1. Let 
us suppose, to begin with, that e = 1. By condition 3' on the basis B% the £ component 
of the spin of the last particle must be ct/v = +1/2. Thus d + (ra, <r) must be at least 1 in 
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this case. Since the value of &n is fixed, the number of states of the basis B\ for which 
d + (n,a) = 5 is now given by (^gl^ 1 )- Therefore 



d(k) 



W) 



B[2d+(n,a)-d(k)] _ n -Bd{k) \^ /"W ~ 1 \ _ 2 B* 

6-1 



^2 q B[2d+(n,a)-d(k)] = q -Bd(k)^ 



5=1 



q B {q B + q' B 



d(k)-l 



(m r = 0, e = 1) 



(4.14) 



If, on the other hand, e = — 1, the only difference with the case e = 1 is that now the £ 
component of the last particle's spin is ctn = —1/2. Obviously, the value of the sum over 
the spins can be obtained from Eq. ()4.14j) by changing the sign of B. We thus have 



q B[2d+(n,a)-d(k)] = g eB fB + g - 



B 



d(k)-l 



(m r = 0) . 



(4.15) 



Inserting Eqs. (|4.13f) and (|4.15|) into the formula (|4.12|) for the partition function we obtain 



Z e (aT) ~ £ 

k&V N 



q B + q- B 



d(k) 



E 



Y, miUi(k) 



mi>-->m r >0 



+ h r ,iq tB [q B + q 



-B 



d(k)-l 



E 



r-l 

Y, miVi{k) 



mi>->m r _i>0 



(4.16) 



It was shown in Ref. [TO] that 

s 

J2 rriiUi(k) 



E 



mi>->m s >0 



-,N, 



3=1 



with 



i=l ^ i=l ' ^ i=\ ' 



(4.17) 



From Eqs. (|4.16|) and (|4.17f) it follows that the partition function Z e of the Hamilto- 
nian (|H.5|) satisfies 



(kl,...,kr)<EV N 



d(k) ( !Ei „ n j 



N i ) VI — g 



N r q eB 
N-+ S kr,l-W 



q a +q- 



. (4.18) 



Substituting ()4.7|) and (j4.18f) into (|4.5|) we finally obtain the following expression for the 
partition function of the Haldane-Shastry spin chain 1)4. 4jl : 



ut) = n 



1 _ Ap+N-\{i+l)) 



E 

(ki,...,k r )eP N 



q B + q~ B 



d(k) 



r_1 „iVi 

.7=1 



1-g 



N r 



r^ q B +q -B 



(4.19) 
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It can be seen that all the denominators 1 — q j appearing in this formula are included as 
factors in the first product. Similarly, the denominator in the last term is always canceled 

by the factor Iq + q \ (indeed, d(k) > 1 for k r = 1). Therefore, as should be the 
case for a finite system, the partition function Z e can be written as a finite sum of terms 
of the form d e q e , where 8e is an eigenvalue of the spin chain Hamiltonian h e and d e its 
corresponding degeneracy. 

The formula (|4.19|) for the partition function Z e reduces to that found in Ref. jlUj in 
the absence of magnetic field. We also note that the spin chain Hamiltonians h + and h_ 
are no longer isospectral when B ^ 0, although their spectra are obviously related by the 
mapping B ^ —B. 
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